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Abstract 

We show that the higher moments of the evolution obtained from the Modified Leading Log- 
arithm Approximation may be regarded as spurious higher order terms in perturbation theory, 
and that neglecting them leads to a good description of the data around and above the peak in 
^ = ln(l/x). Furthermore, we use this study of the moments to show that at high energy the Lim- 
iting Spectrum with Local Parton-Hadron Duality may also be derived from the Modified Leading 
Logarithm Approximation without any non-perturbative assumptions. 
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I. INTRODUCTION 



The perturbative QCD approximation is consistent with a wide range of data. However, 
there are two formal hmitations on the predictive power of this approximation. Firstly, 
perturbative QCD is incomplete in that it does not describe the physics of hadrons entirely. 
Secondly, the perturbation series becomes singular when any of the virtual and real quarks 
and gluons (collectively called partons) in a process has a configuration of energy and mo- 
menta whose "energy scale" E, a quantity whose definition cannot be precisely defined and 
which is integrated over in virtual loops, is as small as Aqcd, the fundamental scale of 
QCD which determines the scale at which non-perturbative effects become important. For 
virtual partons, these two problems are related through the factorization theorem, which 
states that the dominant (leading twist) contribution to a hadronic process at high energy 
is given by a convolution over dimensionless kinematic variables of process dependent quan- 
tities describing those partons with E greater than the factorization scale Q, which is fixed, 
with process independent quantities which describe both asymptotic hadrons and partons 
with E < Q. The former quantities and the Q dependences of the latter quantities are 
calculable in perturbation theory in a given factorization scheme, which defines E, provided 
that Q is sufficiently much larger than Aqcd- Processes with real partons are calculable in 
perturbative QCD provided one performs a physical sum over these particles. 

In the fixed order approach, where the perturbation series is calculated to a finite order, 
the perturbative series for a given partonic process can become divergent in certain regions 
of phase space. Fortunately, the terms that cause the divergences in a given region are 
often calculable to all orders, with the formal sum of all terms of a given class being free of 
divergences in that region and forming the term of a given order in a new series. 

The hadronic and low E partonic components of processes involving the inclusive pro- 
duction of hadrons are contained in fragmentation functions (FF's), which describe the 
probability of transition from a high E parton a to a low E hadron h. In this case it is 
valid to identify E with the transverse momentum of each parton. The cross section for 
a given process is obtained by convoluting the FF's with the partonic cross sections (the 
coefficient functions) over the ratio z of the hadron's longitudinal momentum to that of 
the parton. The FF's at one value of Q can be calculated perturbatively from those at 
another value Qq by convoluting with the evolution matrix, which describes the probabili- 
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ties of transition from a parton at Q to another at Qo, and which is obtained in terms of 
the perturbatively calculable splitting functions by solving the DGLAP equation. However, 
the series for the splitting functions breaks down as 2 — due to terms which behave in 
this limit like (a^/z) In^""^"™ 2, where m = 1, ...,2n — 1 labels the class of terms. (Terms 
which behave like a" are classified as m = 2n.) These logarithms must be resummed before 
the fixed order splitting functions are valid at small z. Since the cross section at hadronic 
momentum fraction x = 2^/1/3, where p is the momentum of the produced hadron h and 
a/s is the centre-of-mass energy, depends on the FF's over the range x < z < 1, such a re- 
summation is required to describe the cross section at small x. Resummation of the leading 
(m = 1) and subleading (m = 2) logarithms, which appear at leading order, is obtained via 
the Modified Leading Logarithmic Approximation (MLLA) (for reviews see js, ^). Since 
the coefficient functions are non-singular as z — > 0, and the quark FF's are proportional to 
the gluon FF in the MLLA, the cross section at low x is then proportional to the MLLA 
evolved gluon FF. 

There is some freedom to choose the MLLA evolution due to the next-to-MLLA error. 
The evolution given by the analytic solution to the MLLA differential equation ^ is well 
behaved for Qq = O(Aqcd); which may imply that the two limitations on perturbation 
theory to describe hadronic physics that were mentioned at the beginning of this section are 
too stringent. The first limitation may be weakened by introducing the Local Parton-Hadron 
Duality (LPHD) hypothesis jsl, which states that the distribution of partons below a certain 
energy scale in sufficiently inclusive processes is similar to the distribution of hadrons, up 
to the number of particles actually produced (the multiplicity). This implies that the initial 
hadronic fragmentation function is proportional to the partonic fragmentation function for 
Qo = O(Aqcd)- The second limitation may be weakened by assuming that partons with 
E = O(Aqcd) may be described by perturbation theory (after resumming with the MLLA). 
Together with the LPHD, the particular choice Qo = ^qcd in the MLLA evolution gives the 
Limiting Spectrum The only free parameters are the initial normalisation and Aqcd, 
which can be fitted to the data. A good fit to the data over the whole range of ^ = ln(l/x) 
can be achieved, however only if the evolution of the normalization is modified. In 6[ an 
additional component not provided by the MLLA was added to the normalization, whereas 
in a different normalization was fitted for each value of ^/s. Otherwise Aqcd obtained in 
this approach is consistent with that of other analyses. 
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In our recent work [8| , we studied the MLL A evolution without using strong assumptions 
such as the LPHD or the vahdity of the Limiting Spectrum, nor modifying the MLLA 
evolution itself. Using an initial scale Qq ^ Aqcd and a parameterized function for the 
initial gluon FF, we achieved a good description of the charged hadron cross section data for 
^ up to and around the peak, and obtained values of Aqcd close to those in the literature. 
Beyond the peak the MLLA evolution turned out not to be sufficient to describe the data. 
The theoretical curves exhibited a second bump after the first peak, not seen in the data, 
which have a characteristic Gaussian shape around the first peak. Such a problem cannot 
be solved by modifying the MLLA normalization. 

It is expected that the fixed order approach with double and single logarithms resummed 
with the MLLA should give a good description of the small to large ^ data, and therefore, 
if fixed order corrections are not included, one has to consider qualitatively what effect the 
fixed order prediction at small ^ has on the MLLA prediction at large ^. In fact, the MLLA 
formally improves the description of the Mellin transform of the cross section for small \uj\, 
where u is defined in Eq. (^, rather than the cross section itself at large ^. Fixed order 
calculations indicate that the large ^ region has rather little dependence on the large |a;| 
region. 

It is the purpose of this paper to consider the features of the data that the MLLA can 
describe, and thereby modify the MLLA evolution in order to improve the description of 
the large ^ behaviour of the spectra without spoiling the description around the peak. We 
start in Section |n] by considering the moments of the cross section, since these quantities 
depend very little on the large \uj\ behaviour of the evolution and the parameterization 
of the initial gluon FF, and can be easily extracted from the data. We then derive our 
approach for improving the large ^ description. In Section IIIIl we compare the predictions 
of this approach with the experimental data at large ^. In Section IIVI we study the effect 
of imposing the limits of the LPHD and Limiting Spectrum on our approach. Finally, in 
Section El we present our conclusions. 

II. EVOLUTION OF MOMENTS IN THE MLLA 

In this Section we outline the features of the MLLA that will be important for the 
derivation of our main result. More details can be found in our previous publication j^, 
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however, to make our formulae here more transparent we will refrain from using the variables 
Y = \ia{Q/Qo) and A = ^(QoMqcd) and write Q and Qo exphcitly. 

The MLLA is believed to describe the energy dependence of cross sections for hadron 
production in the region for which ^ 1 and |ci;| = 0(^/aJ), where uj replaces the variable 
X in the Mellin transform 

POO 

fu = j (i^exp[-u;^]x/(x). (1) 

In this limit the cross section is proportional to the gluon FF D{x,Q) at the conventional 
choice of factorization scale Q = y/s/2. The dependence of the gluon FF on Q in Mellin 
space takes the simple form 



DUQ) = EUas{Q),asiQo))D^iQo) 



(2) 



where D^{Qo) is non-perturbative, while E^^j is determined in terms of the gluon splitting 
function 7^(0:5), 



E^{as{Q),as{Qo)) = exp 



rfln/i 7^(as(/i)) 



(3) 



From the MLLA, the double and single logarithmic contribution to 7a; (as) reads 



H + 



as 
2^ 



u;2 + 47 



o 



uo 



(4) 

where as{Q) is calculated at one loop order and depends on the number of flavours Nf, 
72 = ANcas/{2T:) for A^'^ = 3 colours, a = llNj3 + 2Nf/{3N^) and b = llNj3-2Nf/3. As 
is usual in applications of the MLLA, we choose Nj = 3. Eq. ^ is an expansion of 7^; in 
as/uj keeping ag/uj'^ fixed. The first line is of order agjio^ and is obtained from the Double 
Logarithm Approximation (DLA), while the second is the MLLA correction of O ((as/cu)^). 
The O {[fiLsjiSf^ error is the unknown next-to-MLLA correction. Since Eq. (@]) reduces to a 
finite series in ^Jal at = 0, the MLLA should also be a good approximation in the region 
i/ttl <^ 1 if the next-to-MLLA corrections in this region are of higher order in ^Jo^s. 
One is interested in determining what improvements the MLLA makes to the cross section 
in X space. This quantity can be obtained from the inverse Mellin transform, given by 



xD{x,Q) = / dujexp[uj^]D^{Q), 
2m Jc 



(5) 
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where the contour C may be taken to be a straight hne from uj = ojq — ioo to u = uq + 
ioo, where cuq is real and to the right of all singularities in D^^IQ). Actually, the small 
^ dependence of D{x,Q) is largely determined by at large \u\, which is described by 
the fixed order result. Coincidentally, the MLLA also leads to a good description of the 
small ^ region, since in Eq. (0} becomes negative at large \uj\ like the fixed order result. 
On the other hand, as ^ — 00 the contribution from the \uj\ < 0{y/a^) region becomes 
increasingly relevant, and may eventually become dominant since the evolution of the cross 
section calculated in the fixed order approach falls off rapidly at large |ci;| due to the — In la; | 
behaviour of the anomalous dimensions to all orders. However, 7^ in Eq. only approaches 
a constant, — aas/(27r), at large \uj\, so it cannot be expected that this approach to the 
MLLA gives a good description at large ^. There is no guarantee that there exists some 
suitable choice for the large |a;| behaviour of -Da;(Qo) which can remedy this problem with 
the evolution. 

Therefore we require an evolution which approximates the MLLA well in the \uj\ < 
0{y/Q^) region and which falls off sufficiently fast at large \u!\ such that the small \uj\ region 
gives the dominant contribution to the cross section at large ^. For this purpose we will 
study the MLLA in terms of the moments Kn of the gluon FF, where 



Kn{Q)= (-^) InDUQ) 



(6) 

ui=0 



since the first few moments {n finite) depend very little on the behaviour of the evolution 
at large Ic^l. The K moments completely determine Di^{Q), since Eq. © may be inverted 
using Taylor's Theorem to give, formally. 
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\nD^ = J2^^I^. (7) 
^ — ^ nl 

n=0 

Note that the K moments may be expressed in terms of the normalized ^ moments, (■C"), 
which from Eq. (Q) can be calculated using 



(8) 



0^=0 



From Eqs. and (jH)), the moments of the gluon FF evolve as 

KniQ) = KniQo) + AKniasiQ), a.(Qo)), (9) 

where AKn is defined to be the nth K moment of the evolution E^^, and so the K moments 
have the additional advantages that their MLLA evolution is independent of D^{Qq) and 



that they each evolve independently of one another. This definition of Ai^„ together with 
Eq. Q implies 



(10) 

uj=0 



AKn{as{Q),as{Qo)) = J rfln/x |^-— J 'y^{as{fi)) 
Explicitly, Eqs. ^ and (fTUI) for n > 1 give 
AK„(a,(Q),a,(Qo)) = a;^{Q) (cf + C«a|(g) + O (a,)) - {a,(g) ^ a.(go)}, (11) 



where the O (a^) error refers to unknown next-to-MLLA corrections, while the Cn'^^ are 
completely determined and are presented in ^ for the first few values of n. In fact, for 
n > 3 and odd, Cn^ = 0. Eq. also applies for n = 0, but AKo also contains a term 
proportional to Ina^. 

For small \u!\, E^{as{Q),(ys{Qo)) may be approximated by 

^ f (n\ (n\\ A/s:„(a,(g), a,(Qo))(-^)'' 
ln£;^(a,(Q),a,(Qo)) = > j (12) 

n=0 

with M finite. This corresponds to evolving the moments Kn for n < M exactly as in the 
full unexpanded case, but fixing the remaining K moments to be equal to Kn{Qo)- Since 
the coefficient Cn'' in Eq. (fTT|) for n even is positive when n/2 is odd and negative when 
ri/2 is even, then in the region where the imaginary part of uo is large, if M > 2 and even, 
Q and Qq are large and Qq < Q, we obtain the fast decrease 



exp 



Ml 



(13) 



If M is odd, the Mth term in Eq. p2|l just produces an oscillation, in which case the 
replacement M — >■ M — 1 in Eq. (jl3|) must be made. 

It remains to be found whether Eq. (fT^ agrees well with the approach of Eqs. Q to (jlj) 
in the whole region \uj\ < 0{y^). In the extreme case \uj\ = 0{^ya^), Eq. (fTTj) implies 

Lj^AKr, = O (a;'A , (14) 



so that all terms in the series in Eq. (jl2|) become of similar magnitude. Such a series may 
still converge, or oscillate with an average value equal to the unexpanded result. In any 
case, we see that the accuracy of Eq. (fT^ to reproduce the MLLA contribution to the cross 
section cannot be reliably determined in Mellin space. Therefore we will try to determine 



what the suppression of higher moments means in x space. For this purpose, it will be 
convenient to work with the moments A^, S,-, cr^ and for n = 3, oo, defined by 

, 3, 



a-" 



K4 



((e-e)') 



3, Kji 



(15) 



{n > 5). 



Note that K3 is often written as s and K4 as k. From Eqs. (jHI) and (jH)), these moments are 
related to the K moments by 



Ko = lnN, K, = ^, K2 = (y\ ir„ = a"«:„ (n > 3). 



(16) 



To obtain an expression for some function D{x) in terms of the we first make the 
replacement y = iua in Eq. (0), which yields 



xD(x) 



N 



exp 



(17) 



where 5 = (^ — 0/^ ^ind the real quantity R is given by 



foo 



dyexp 



.n=3 



exp 



iy5 



y_ 

2 



21 



(18) 



Note that R is equal to unity when all the k„ vanish. We therefore see that a function 
depends on its moments Kn for n > 3 only through the ratios k„ = Kn/cr"'. To calculate Eq. 
(II 7|) to a given accuracy when the k„ are small, we expand the K„-dependent exponential in 
Eq. (jl8j) in powers of the k„ up to the required accuracy and perform the integral for each 
term. R is then rewritten as an exponential of the form 



R = exp 



i=0 



(19) 



where each Ai vanishes when all the k„ vanish. Now suppose the k„ are sufficiently small 
such that it is valid to expand the Ai in the k„ for n < M up to some finite order, and 
neglect the moments «:„ for n > M. In this case In i? as a series in 6 will terminate at some 
finite order, hence the argument of the exponential in Eq. (fT^ is not an expansion in S, 
since the approximation is valid even if 6 is of 0(1). For example, including the complete 
contribution from all terms of 0(s), 0{k), O(s^), 0{k'^), 0{sk), 0{k^) and 0{kq) in Ini? 
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gives 



xD(x) 



N 



exp 



-fc --sb--{2^ k) 6^ + -s6^ + ^kS^ 
8 2 4 ^ ^ 6 24 



5 2 1,2 1 /I 2 ,\ /I 2 1,2 1 \ .2 

H fc^ Kg + -K5 sk ] 6 + { -s^ k^ H Ke 6^ 

24 12 48 VS 3 J \2 3 16 ' 

+ f.l^^ + 53 ^ /^_i^2 1^2 _ 1 \ ^4 

V 12 12 / V 8 48 48 / 



(20) 



+ 



120 12 / V 72 720 ' 



We now return to hadron production cross sections in the MLLA. In x space, Eq. (j21) 
becomes 

xD{x,Q)= —-E(-,as{Q),as{Qo))zD{z,Qo), (21) 

J X Z Z \ Z / 

where E{z,as{Q),Ois{QQ)) is the inverse Melhn transform of Ei^{as{Q), Ois{Qo)). From Eq. 
(fTHIl . the K moments of zE{z ^ as{Q) , asiQa)) , /tf (as(Q), ^^(Qo)), obey (omitting arguments 
for brevity) 

We may use Eq. (fTT|l to expand (^^(Q), a;s((5o)) as a series in as{Q) keeping asiQo) fixed, 
giving 

o^_0 I- /I -,-1 

(23) 



«:f(a,(Q),a,(Qo)) oca.* (Q) [l+0(^a|(g) 
Therefore we may treat the as small, in which case E takes the form of Eq. (fTTjl . 



2:^(2;, as((5),as((5o)) 



(T^V2^ 



exp 



^"('5"(^),{<}), 



(24) 



where 5^{z) = (ln(l/2;) — ^ )/o"^- Eq. and the discussion after Eq. (fT^ show that 
when Ini?-^ is expanded in as{Q) while 6^{z) is fixed, the higher moments serve only to 
contribute spurious higher order terms to this series. Such terms give a large theoretical 
error, since they contribute unstable "noise" , as can be seen by performing the evolution on 
a smooth function. Discarding these terms gives an evolution in which the higher moments 
of the gluon FF are fixed with respect to Q, in other words Eq. (fT^ . 
Similarly, the k moments of D{x,Q) obey 

KniQo) + AKn{as{Q),as{Qo)) 



Kn{Q) 



[K^iQo) + AK2{as{Q),a,{Qo))]- 



(25) 
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TABLE I: Simultaneous fit of Aqcd and the moments InA^, ^, o"^, and at Qq = 14/2 GeV 
to the same moments of the data from Xdf ~ 



IniV 






K4 


Aqcd (MeV) 


2.00 


2.09 


0.40 0.05 


-0.19 


411 ± 36 



At sufficiently large Q, AKn{as{Q), as{Qo)) > Kn{Qo), so that, from Eqs. (j22I) and (j23), 
Kn{Q) may be approximated by K^{as{Q),(ys{Qo)), and hence may be treated as small. 
Thus if Qo is sufficiently large, xD{x,Qo) may be parameterized as a distorted Gaussian 
around the average value of ^. 

III. FITTING TO DATA 

Since the evolution of the first few K moments is independent of the approach used (e.g. 
the approach of Eqs. Q to or the approach of Eq. p2j) ) and the parameterization of 
the non-perturbative input, the ability of the MLLA to describe the data in principle can be 
determined by comparison to the moments of the data, assuming that fixed order corrections 
can be neglected. We perform a single fit of Aqcd and the initial Kn with Qq = 14/2 GeV 
to the Kn of the experimental data at various y/s by setting Q = y/s/2 and evolving in the 
MLLA as in Eq. (|TO|l . We use the moments calculated in p], which are presented in the 
form A^, ^, cr^, s and k, from which we extract the first 5 K moments, and their errors are 
obtained by differentiating and adding in quadrature. The results are shown in Table IJ and 
Fig. H The error on Aqcd in Table H] and in all other tables in this paper is obtained by 
inverting the correlation matrix, which is identified with the matrix of second derivatives. 
In Fig. ^ we see that the first 3 moments are fitted very well, however there is a marginal 
disagreement of K3 and K4 with the data which may result from either an inability of the 
MLLA to describe these higher K moments or from the inaccuracy involved in obtaining 
these moments from the experimental data, in which case the true experimental error would 
be larger than that shown. We also perform a fit directly to the basis of moments used in , 
and found no significant change in the values of the initial parameters, and the theoretical 
curves for s and k deviated seriously from the data below y/s ~ 25 GeV. 

These results strongly suggest that there exists some approach to applying the MLLA and 
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FIG. 1: Simultaneous fit of tlie initial moments and Aqcd to all of the experimental data moments 



calculated in 



, by evolving the moments in the MLLA. 



some parameterization of the non-perturbative components that gives a good description 
of the data over a large range in ^. In the following three subsections we perform fits 
directly to the data points at different ^ using various approaches. In Subsection 1111 Al 
we evolve the moments, place them in a distorted Gaussian and compare with the data 
at different ^. In Subsections IIII Bl and IIII TJ\ we evolve in Mellin space using Eq. (fT^ . 
with two different parameterizations of the non-perturbative input: In Subsection IIII Bl we 
parameterize the initial distribution xD{x,Qq) such that its higher moments are exactly 
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zero but the remaining moments are left as free parameters, while in Subsection IIII CI we 
parameterize the initial distribution as a distorted Gaussian (the first two lines of Eq. (HOJ) 
so that the higher moments are small. In order to avoid the small ^ region, where fixed 
order effects are important and where the data have a high accuracy, we follow |7|] and use 
only those data for which 

^ > 0.75 + 0.33 In(v^), (26) 
but impose no upper limit in ^ on the data. 



A. Evolution of moments as distorted Gaussian parameters 



While the data tend to follow the shape of a distorted Gaussian well, the distorted 
Gaussian evolved with the approach of Eqs. (0) to (jH) does not — in jQ], the theoretical 
curves exhibited two bumps. Therefore we constrain the evolved xD{x,Q) to follow a 
distorted Gaussian, 



xD{x, Q) 



N' 



aW27T 



exp 



1 1 
-k' - -s'6' 
8 2 



i (2 + k') 5" + + j^k'5" 



(27) 



with 6' = {C, — ^)/o-', where the parameters A^', ^ , a', k' and s' depend on Q. Since these 
parameters are approximately equal to the corresponding unprimed quantities defined in 
Eq. (fT^. we choose their Q dependences to be the same, i.e. that obtained from Eqs. (fTU)) 
and (fT^ . By comparing xD[x,Q) calculated in this way with the data at different ^ and 
y/s, we fit the initial K'^{Qq) = K'^ for n < 5, as well as Aqcd, with the choice Qq = 14/2 
GeV. A similar approach was applied in Q, however the Limiting Spectrum formulae for 
the moments were used, and Aqcd, the peak position of the data and the normalisation of 



the data^for each y/s were fitted. We use TASSO data at y/s =14, 22, 35 _ and 44 GeV 1^|, 

3, 



TPC Q and MARK II LJ data at 29 GeV, TOPAZ data at 58 GeV 
DELPHI Q, L3 161, OPAL [17| and SLD [18l data at 91 GeV, ALEPH 
20| data at 133 GeV, DELPHI data at 161 GeV Q, OPAL data at 172, 183 and 189 GeV 
22 1 and OPAL data at 202 GeV The results are shown in Table ITTj and some of the data 
with the corresponding fitted theoretical curves are shown in Fig. |21 (In all plots in this 
paper, each curve is shifted up from the curve below by 0.8 for clarity.) The fit is excellent 
around the peak and above, but poor below the peak region, where fixed order effects are 
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ALEPH 
131 and OPAL 
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TABLE II: Simultaneous fit of Aqcd and tlie moments IniV', a'^, and K'^ at Qq = 14/2 GeV 
to the data in x space by evolving them in the MLLA and placing them in a distorted Gaussian. 
xIf = 2-4. 



InN' 


r 




Aqcd (MeV) 


2.41 


2.07 


0.91 -0.58 


-0.66 59 ± 4 



important. In particular, the curves show that the MLLA alone predicts the evolution of 
the normalization very well at large C,, in contrast to other analyses. This procedure has a 
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FIG. 2: Global fit by using a distorted Gaussian in which the moments are evolved in the MLLA. 
The lower limits of the data used, given by Eq. ()26j) . are indicated by vertical dotted lines. Each 
curve is shifted up by 0.8 for clarity. 

number of features that differ from that used in . The evolution used there coincident ally 
shares similar properties at large \uj\ with the fixed order result, and hence a good fit was 
obtained with the data below the peak while a large disagreement was found above the peak. 
However, no properties of the evolution used for the fit of Table HTl are shared with the fixed 
order result below the peak, and only the first 5 moments are evolved while the remaining 
moments are held fixed (within the accuracy of the distorted Gaussian approximation). To 
ensure that the deviation below the peak was not due to the lack of data being fitted there, a 
fit using all data for which ^ < ln(-yi/20.5 GeV) was performed, giving a bad fit everywhere. 
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TABLE III: Simultaneous fit of N, ^, cj^, s, k, K5, kq and Aqcd to TASSO data at 14 GeV and 
OPAL data at 91 and 202 GeV, with only the first 3 moments evolved according to the MLLA. 
X^F = 1-14. 

N I a"^ s k K5 Aqcd (MeV) 

10.90 2.40 1.05 0.05 0.54 0.09 0.19 66 ± 10 



TABLE IV: As in Table lnT| but with only the first 5 moments evolved according to the MLLA. 
Xdf = 0-60. 

N I s k K6 Aqcd (MeV) 

10.58 2.15 0.79 -0.79 -0.40 -0.73 -2.44 81 ± 11 



In particular since the parameters try to fit to the accurate data below the peak, a large 
initial {K^l was obtained. 

B. Mellin space parameterization and evolution 

We now constrain the initial distribution such that Kn = exactly for n > 7. From Eq. 
((7j), this means that the parameterization of the initial distribution in Mellin space must 
take the form 

,„„„,g„).X:™^. (28) 

n=0 

We perform fits to TASSO data at 14 GeV and OPAL data at 91 and 202 GeV, cut according 
to Eq. and evolve according to Eq. ^ with Qo = 14/2 GeV. Taking M = 2, 4 and 
6, we obtain the results shown in Tables ITTTl II VI and IVI respectivelv. Note that using M = 4 
gives the lowest Xdf- For the case M = 6, we get a large result for kq, since although 
AKQ^asiQ), (^s{Qo)) becomes positive for Q ^ 00, for the data used AKq is negative. This 
is due to the MLLA term being larger than the DLA term, so it may be the case that 
corrections beyond the MLLA are required for this quantity, or some other approach. At 
any rate, from the end of SectionlHl AKq contributes spurious higher order terms and should 
be neglected, and this is confirmed in Table IVl 
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TABLE V: As in Table 11111 but with only the first 7 moments evolved according to the MLLA. 
X^F = 0.77. 



N 


I s 


k 


K5 


KG Aqcd (MeV) 


10.06 


1.80 0.474 -4.86 


2.00 


-1.17 


113.57 96 ± 1 



TABLE VL Simultaneous fit of N', C , cJ'^ s', k' and Aqcd to TASSO data at 14 GeV and OPAL 
data at 91 and 202 GeV, with only the first 3 moments evolved according to the MLLA. Xdf ~ 1-67. 



N' 


r 




s' k' 


Aqcd (MeV) 


10.51 


2.22 


1.17 


-0.71 1.18 


86 ± 1 



C. Distorted Gaussian w^ith Mellin space evolution 

In global analyses, the initial FF's are parameterized in x, and the parameters are fitted 
to data by evolving the FF's in Mellin space in the fixed order approach. For this reason, 
in ^ the initial gluon FF was parameterized as a distorted Gaussian and evolved using the 
approach of Eqs. to (jH). We repeat this approach here, but instead we will apply MLLA 
evolution in the form of Eq. (fT^ . First we repeat the fits of Subsection IIII Bl The results 
for the case M = 2 and M = 4 are shown in Tables IVII and I VI II respectively. Note again 
that using M = 4 gives the best fit. The resulting curves for the M = 4 case are shown 
in Fig. El The cases M = 6 cannot be tested since, as we found in Subsection IIII Bl AfCg 
is negative for the data used so that the integral for the inverse Mellin transform does not 
converge. Fits in which data for all is used generally give a bad fit everywhere except at 
values of ^ beyond the peak, which is due in part to the high accuracy of the data at small 
^. As anticipated from the fit of Table m excellent agreement is found around and above 
the peak region. 

We now perform a fit to all the data, using M = 4 in the evolution. The results are 
shown in Table IVinl and Fig. 01 and are the main results of this paper. 

In all our approaches, reasonably consistent values of Aqcd were obtained of around 100 
MeV. The small ^ region was generally not well described, however this region is outside 
the scope of the MLLA and requires fixed order corrections. It is generally found that the 
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TABLE VII: As in Table IV1[ but with only the first 5 moments evolved according to the MLLA. 
X^F = 0.65. 



N' 


a'2 


s' k' 


Aqcd (MeV) 


10.75 2.12 


0.92 


-0.65 -0.22 


77 ± 10 



15 




' ' ' ' ' ' ' ' ' ' 1 

1 2 3 4 5 6 7 



FIG. 3: Fit of the distorted Gaussian parameters and Aqcd to TASSO data at 14 GeV and OPAL 
data at 91 and 202 GeV, with only the first 5 moments evolved according to the MLLA. 

best fit is obtained with M = 4, although with other values of M we also obtained good fits 
around and above the peak. 



IV. THE LPHD AND LIMITING SPECTRUM LIMITS 

From Eq. (fTT|) . AKn{as{Q), as{Qoj) ^ oo as Q — > oo, so that for sufficiently large Q the 
initial Kn{Qo) in Eq- & may be neglected. However, such an approximation should not be 

TABLE VIII: Global fit with only the first 5 moments evolved according to the MLLA. Xy)F ~ 



N' 


r 




s' k' Aqcd (MeV) 


10.39 


1.90 


1.10 


-1.45 -0.09 114 ± 6 
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FIG. 4: Global fit with only the first 5 moments evolved according to the MLLA. 

made for n = 0, since the cross section is very sensitive to the initial InA^. Therefore, data 
at sufficiently large Q should be reasonably well described with a starting distribution of 
the form xD{x,Qq) = N6{1 — x). Thus, from a purely perturbative analysis, we see that 
the LPHD arises because the perturbative components form the dominant contribution to 
the cross section. In this sense, the LPHD follows from the MLLA rather than being an 
additional assumption. 

For sufficiently large Q, any term of 0{a~'^{Qo)) with n positive can be neglected relative 
to a term of 0{a~'"'{Q)), so that from Eq. (fTTj) . we may replace AKn{as{Q),C(s{Qo)) with 
AKn{as{Q), oo). This can be artificially achieved by setting Qq = Aqcd, since then as{Qo) = 
oo as a consequence of the perturbative approximation. Thus the Limiting Spectrum also 
follows from the MLLA. However, it is important to note that the AKn{<ys{Q),oo) for n 
sufficiently small would not be finite if corrections of next-to-MLLA or higher were included. 

To summarize, a simple analysis of the MLLA shows that, at sufficiently large Q, the 
assumptions of the LPHD and the Limiting Spectrum will appear to be correct, since the 
cross section may be calculated with Eq. Q approximated for n > 1 by 

i^„(Q)^Air„(a,(g),oo). (29) 

To properly test the physical assumptions which imply the LPHD and the Limiting Spectrum 
requires a comparison with data in the region Q = O(Aqcd)- However, the only essential 
difference between our procedure and that of the LPHD with the Limiting Spectrum is 
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that we do not allow for any assumptions on the size of the initial ^ and a^, we only 
assume that the initial k„ are small. We now study what the effects of these additional 
constraints are, by using Eq. ()29|1 for n > 1. From Eq. 0, the normalization N{Q) obeys 
InA^(Q) = lnN + AKo{as{Q),oo), where IniV = \nN{Qo)-AKo{as{Qo),oo) is independent 
of Qq. Thus our predictions will be exactly independent of Qo, and the theory contains just 
two free parameters, and Aqcd, which we fit to all available data. We find that the data 
above the cut is best described when the cut is chosen as in Eq. but with = +1 
added to the right hand side. Fitting a different normalization for each y/s of the data does 
not improve the fit significantly. For M = 4 in the evolution, we obtain = 7.68 ± 0.02, 
Aqcd = 446 ± 2 and Xdf = 9.3, while for M = 2 we obtain iV = 6.64 ± 0.02, Aqcd = 292 ± 1 
GeV and Xdf ~ 3.74. The resulting plots are shown in Fig. El for the M = 2 case. (The 
M = 4 fit gives similar curves.) With these results we are able to calculate the moments 
at Q = 14/2 GeV, and find reasonable agreement with our previous results for A^, ^ and 
(e.g. that of Table rvTTT|) . The evolution of these quantities follows the data, as is to be 
expected since we are using the same evolution as we used in the fit of, e.g.. Table IVlTTl 
However, it is clear that the prediction at and below the peak is significantly dependent 
on the initial values of ^ and cr. Comparing with other analyses which use the LPHD with 
the Limiting Spectrum, we see that the data above the peak are well described, with no 
modifications to the normalization, while the data below the peak are not. 

15 I ' , ' , ' , ' , ' , ' , ' 1 




1 2 3 4 5 6 7 

FIG. 5: Global fit to data with the LPHD and Limiting Spectrum, and with only the first 3 
moments evolved. 
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V. CONCLUSIONS 



In j8| , it was shown that a naive apphcation of the MLLA without additional assumptions 
gives a good description of data around the peak region, but not beyond. Since the evolution 
used approached a constant at large |ci;|, while the fixed order approach implies that the 
evolution falls to zero, the contribution to the cross section at large ^ from the cross section 
at small Ic^l was underestimated. This qualitative feature of the fixed order approach can 
be taken into account by using the Limiting Spectrum, which decreases fast at large Ic^l. 
In this paper, we found that this feature can also be taken into acount by suppressing the 
evolution of the higher moments. This does not affect the approximation since, by studying 
the MLLA in x space, one can see that it is formally justified to neglect their evolution 
in the perturbative approximation. In addition, at sufficiently large Q, one finds from the 
MLLA that the spectrum acquires a distorted Gaussian shape. Fixing the higher moments 
and using a distorted Gaussian for the initial distribution at Qq = 14/2 GeV resulted in 
a good description of all data for which > 14 GeV from just below the peak to the 
largest value of ^, and one obtains Aqcd ~ 100 GeV. In obtaining these results, we used 
a non-perturbative input that was determined empirically at a low scale, rather than from 
physical arguments such as those of the LPHD. 

Furthermore, we showed that the cross section approaches that of the Limiting Spec- 
trum at sufficiently large Q. We stress that this follows from the MLLA, without additional 
hypotheses. In practice, imposing the Limiting Spectrum limit on the evolution with sup- 
pressed moments and imposing the LPHD on the initial distribution gives a poor description 
up to the peak. However, a good description is obtained beyond the peak without requir- 
ing any modification to the normalization, and the best fit, where M = 2, resulted in 
Aqcd = 292 GeV, which is consistent with other analyses. 

Finally, the region of the data used in global fits may be extended to lower values of x 
(large ^) by incorporating the MLLA into the fixed order calculations. This would allow 
for a better determination of the FF's, particularly at small momentum fractions, as well as 
more constraints on Aqcd- 
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